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Abstract 

We present measurements of the transverse momentum distribution of W 
and Z bosons produced in pp colhsions at \/s = 1.8 TeV. The data were 
cohected with the D0 detector at Fermilab during 1994-1996 Tevatron run. 
The results are in good agreement with theoretical predictions based on the 
perturbative QCD and soft gluon resummation combined calculation over the 
entire measured pT range {pT = — 200 GeV/c). 
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I. INTRODUCTION 



At the Fermilab Tevatron, W and Z bosons are produced in high energy pp coUisions. The 
study of the production of W and Z bosons provides an avenue to explore QCD, the theory of 
strong interactions. The benefits of using intermediate vector bosons to study perturbative 
QCD are large momentum transfer, distinctive event signatures, low backgrounds, and a 
well understood electroweak vertex. In this paper we present measurements of the W and Z 
transverse momentum distribution based on data taken by the D0 collider detector during 
the 1994-1996 Tevatron running period. 

In the parton model at lowest order, W and Z intermediate vector bosons are produced in 
head-on collisions of qq constituents of the proton and antiproton, and have little transverse 
momentum {pt « M^, Mz)- Consequently, the fact that observed bosons have large 
transverse momentum {pj) is attributed to the production of one or more gluons or quarks 
along with the bosons. At high transverse momentum {px > 20 GeV/ c), the cross section is 
dominated by the radiation of a single parton with large transverse momentum. Perturbative 
QCD is therefore expected to be reliable in this regime |l|]. At low transverse momentum 
{Pt < 10 GeV/ c), multiple soft gluon emission is expected to dominate the cross section. A 
soft gluon resummation technique is therefore used to make QCD predictions. Neither 
the resummed nor the fixed order calculation describes the distribution for all values of pt- 
Conventionally, one switches from the resummed calculation to the fixed-order calculation 
aX Pt ^ Q [iiil- Thus, a measurement of the transverse momentum distribution may be 
used to check the soft gluon resummation calculations in the low px range, and to test the 
perturbative QCD calculations at high px- 

II. INTRODUCTION TO THEORY 

In standard fixed order perturbative QCD, the partonic cross section is calculated by 
expanding in terms of the strong coupling constant (a^) [|l|, where each power of as cor- 
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responds to the radiation of a single gluon or quark into the final state. This procedure 
works well in the high p-r region, where p"^ ~ Q^. However, as pt 0, this fixed order 
calculation of the cross section diverges due to the presence of correction terms that go as 
log^ {Q"^ / p^) . Physically, this failure of fixed order perturbation theory at low pt is due to 
soft gluon radiation from the initial partons that is not properly accounted for in the stan- 
dard expansion. This difficulty in performing the calculation can be remedied by reordering 
the perturbation series through a technique called resummation, where the cross section is 
calculated in terms of the large logarithms, rather than strictly in terms of powers of 
The resulting calculation is an all orders calculation, i.e., each piece in the new sum contains 
terms to all orders in a^. However, the largest terms dropped in the latest calculation are 
0{a^), so it is considered to be accurate to O(a^). 

In final form, the calculation is carried out via a Fourier transform in impact parameter 
space (6-space), with the following relation describing the differential cross section [2-5]: 



d 



'2, 



d% e'f"^-' W{b, Q) + Y{b, Q) (1) 



dp'^dy 

where W{b, Q) contains the results of resumming the perturbative series, and Y{b, Q) adds 
back to the calculation the pieces that are perturbative in a^, but not singular at = [Q. 

Although the resummation technique extends the applicability of perturbative QCD to 
lower values of px? a more fundamental barrier is encountered when pt approaches Aqqd- In 
this region non-perturbative aspects of the strong force dominate the production of vector 
bosons and in general perturbative QCD is expected to fail. This implies that the resummed 
calculation becomes undefined above some value of 6 = bmax- In order to extend the calcula- 
tion to the low Pt region a parameterization which accounts for the non-perturbative effects 
is introduced. This extension is accomplished by cutting off the integral in Eq. |1] at some 



value bmax and replacing W{b, Q) with W{b^, Q)e~^'^p^^''^\ where 6* = b/ ^Jl + b/b^ax- This 
effectively cuts off the contribution of W{b, Q) near bmax, leaving the cross section dominated 
by the function being introduced SNp{b,Q), known as the non-perturbative Sudakov form 
factor. Snp has the generic renormalization group invariant form [||: 
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SNp{b, Q) = h{b, xa) + h{b, xb) + h^ib) In (g) (2) 

where xa and xb are the momentum fractions of the incoming partons, h is Fourier conjugate 
to the transverse momentum (impact parameter), Qo is an arbitrary momentum scale and 
hi{b,x) and h2{b) are phenomenological functions to be determined from experiment 
We used the Sudakov factor functional form from Ladinsky and Yuan ||^: 

Sji^p{b,Q) = g,b^ + g2bHn[^) + g,b\n{lOOxAXB) (3) 

The values of Qi were determined by Ladinsky and Yuan by fitting to low energy Drell- 
Yan data and a small sample of Z ^ ee data from 1994-96 run at CDF 0, yielding gi = 
O.lltn^ GeV2, g2 = OAtoi GeV^ and g^ = -1.5tl]i GeV-\ where b^ax = 0.5 GeV~^ and 
Qo = 1.6 GeV, and the CTEQ2M pdfs were used. 



III. MEASUREMENT OF THE DIFFERENTIAL CROSS SECTIONS 

In this paper we present measurements of the pj- spectra of W and Z bosons produced 
in pp collisions at = 1.8 TeV with the D0 detector at Fermilab. The transverse 
momentum spectra of W and Z bosons have been measured previously by the UAl [|^, 
UA2 [|], CDF and D0 |jlO| collaborations, but with smaller data samples than the ones 



reported on here. 

The W and Z samples have been selected from data taken during the 1994-96 run of the 
Tevatron, and correspond to an integrated luminosity of ~ 80 pb^^ for W and ^110 pb^^ for 
Z. The measurements of the W and Z boson p^ spectra used the decay modes W ^ eu and 
Z e~^e~ . Electrons were detected in hermetic, uranium liquid-argon calorimeters with an 



energy resolution of about 15% / y E(GeV) . The calorimeters have a transverse granularity 
of Ai] X A(f) = 0.1 X 0.1, where rj is the pseudorapidity and (p is the azimuthal angle. 
Electrons were accepted in the regions \ri\ < 1.1 (central) and 1.5 < |?7| < 2.6 (forward). 
In reconstructing the px of the W boson, we assume that the transverse momentum of the 
neutrino is given by the calorimetric measurement of the missing transverse energy (^j, ) in 



the event. Electrons from W and Z boson decays tend to be isolated. Thus, we required 
the cut 

Et,,{OA) - Eem{0.2) 

where £"^4(0.4) is the energy within AR < 0.4 of the cluster centroid {AR — Ar]'^ + A0^) 
and Eem{0-'2) is the energy in the EM calorimeter within AR < 0.2. 

At trigger level, the Z analysis required two electrons, one with transverse energy (Et) 
greater than 20 GeV and the second with Et greater than 16 GeV. Off-hne for a "loose" 
electron we required that Et > 25 GeV and that the transverse and longitudinal shower 
shapes be consistent with those expected for an electron (based on test beam measurements) . 
For a "tight" electron we also required a good match between a reconstructed track in the 
drift chamber system and the shower position in the calorimeter. For the Z boson sample 
we required one electron to be "tight" and the other to be either "tight" or "loose" ; at least 
one of the two electrons had to be in the central region. To be acceptable candidates for Z 
production, both decay electrons are required to be isolated. The dielectron invariant mass 
was required to lie in the range 75 — 105 GeV/c^. These selection cuts were passed by 6407 
Z candidates. 

In the case of the W boson analysis, a single electron with Et greater than 20 GeV 
was required at trigger level. Off-line, a tighter requirement on the electron quality was 
introduced to reduce the background level from QCD events, especially at high transverse 
momentum. Electron identification was based on a likelihood technique. Candidates were 
first identified by finding isolated clusters of energy in the EM calorimeter with a matching 
track in the central detector. We then cut on a likehhood constructed from the following four 
variables: the from a covariance matrix which measures the consistency of the calorimeter 
cluster shape with that of an electron shower; the electromagnetic energy fraction (defined as 
the ratio of the portion of the energy of the cluster found in the EM calorimeter to its total 
energy); a measure of the consistency between the track position and the cluster centroid; 
and the ionization dE/dx along the track. For the W boson sample we required one isolated 
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electron in the central region and Ipj. > 25 GeV. The event is rejected if there is a second 
electron and the dielectron invariant mass lies in the range 75 — 105 GeV/c^. A total of 
41173 W candidates passed these cuts in the central region. 

The absolute normalization of the trigger and selection efficiencies were determined using 
Z — i> e~^e^ D0 collider events in which one of the electrons satisfied the trigger and selection 
criteria. The second electron then provided an unbiased sample with which to measure 
the efficiencies. The variation in the electron selection efficiency as a function of pt has 
been determined using a signal Monte Carlo sample from HERWIG, smeared for detector 
resolutions and overlayed with events taken in random pp collisions (zero bias). A parametric 
Monte Carlo program |Tl[] was used to simulate the D0 detector response and calculate the 
kinematic and geometric acceptance as a function of px- The detector resolutions used 
in the Monte Carlo were determined from data, and were parameterized as a function of 
energy and angle. The relative response of the hadronic and EM calorimeters was studied 
using the transverse momentum of the Z boson as measured by the pt of the two electrons 
compared to the hadronic recoil system in the Z event. This parameterized representation 
of the D0 detector was used to smear the predictions by detector effects and compare it to 
our measured results. 
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FIGURES 
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FIG. 1. Left: The background fraction to the W ^ ev signal from QCD, Z — > ee , and top 
production as a function of the px of the W . Right: The background fraction from QCD processes 
to the Z — > ee signal as a function of the pT of the Z. 

The major source of background for both the W and the Z sample is QCD multi-jet and 
photon-jet events: the amount of background in the samples and its shape as a function of 
Pt was obtained directly from D0 data. 

The total amount of QCD background in the inclusive W sample is ~ 1%. The nor- 
malized QCD background was subtracted bin-by-bin from the W boson candidate sample 
transverse momentum spectrum. Additional corrections were made to account for top quark 
background events (0.1%) and for Z — > e^e~ events (0.8%), where one of the electrons was 
lost or not identified. Since was measured from the recoiling hadrons, the events origi- 
nating from W ^ TV (where r — >■ evv) contributed properly to the differential distribution; 
this source of background therefore was included in the Monte Carlo simulation of the 
distribution. 

The total background for the Z sample is ~ 2% for events in which both electrons are in 
the central region and ~ 7% when one electron is in the central region and the other is in the 
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forward region. Backgrounds to Z — > ee production from Z tt, top quark and diboson 
production have been estimated from Monte Carlo samples and are negligible. Figure |1| 
shows the background fractions as a function of px for both the W and Z samples. 




FIG. 2. Left: D0 result (solid points) with statistical uncertainty. The theoretical cal- 
culation by Arnold-Kauffman has been smeared for detector resolutions and is shown as two 
lines corresponding to the ztlo" variations of the uncertainties in the detector modeling. Data and 
theory are independently area normalized to unity. The fractional systematic uncertainty on the 
data is shown as a band in the lower portion of the plot. Right: The ratio (Data-Theory)/Theory 
shown as a function of with its statistical uncertainty as error bars. The total systematic un- 
certainties, shown as a band, were obtained by adding in quadrature the contributions from data 
(background and efficiency) and from the detector modeling. 

The result for the W px distribution, shown in Figure |^, is compared to the theoretical 
calculation by Arnold and KaufFman 0, smeared by detector resolutions. The W data 
shows good agreement with this combined QCD perturbative and resummation calculation 
over the whole range of px- In the case of the Z, we correct the measured cross section 
for the effects of detector smearing. Figure ^ compares the final, smearing-corrected Z 
Pt distribution to the calculation by Ladinsky and Yuan 0. In addition. Figure]^ compares 
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the measurement to the fixed-order perturbative theory |T|. In comparing to the NLO 
calculation, we observed strong disagreement at low p^, as expected due to the divergence 
of the calculation at pt = 0, and a significant enhancement of the cross section from soft 
gluon emission. 
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FIG. 3. Left: D0 Z differential cross section (circles) as a function of pT compared to 
the calculation by Ladinsky-Yuan P|. Theory is normalized to the data. Right: The ratio 
(Data-Theory) /Theory function of pT of the Z. 
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FIG. 4. Left: D0 Z pT smearing-corrected data (solid points) with total uncertainty shown 
compared to the fixed-order (NLO) perturbative prediction The data is normalized to the D0 



measured inclusive Z production cross section |12|; the theory is normalized to its own prediction 
Right: Fractional difference between the Z data and the fixed-order calculation as a function of pT- 
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IV. CONCLUSIONS 



Using data taken with the D0 detector during the 1994-1996 Tevatron collider run, we 
have presented measurements of the W and Z transverse momentum distributions that agree 
well with the combined QCD perturbative and resummation calculations. 
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